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Abstract Are solids intrinsically different from liquids? Must a finite stress be applied in 
order to induce flow? Or, instead, do all solids only look rigid on some finite timescales and 
eventually flow if an infinitesimal shear stress is applied? Surprisingly, these simple ques- 
tions are a matter of debate and definite answers are still lacking. Here we show that solidity 
is only a time-scale dependent notion: equilibrium states of matter that break spontaneously 
translation invariance, e.g. crystals, flow if even an infinitesimal stress is applied. However, 
they do so in a way inherently different from ordinary liquids since their viscosity diverges 
for vanishing shear stress with an essential singularity. We find an ultra-slow decrease of the 
shear stress as a function of the shear rate, which explains the apparent yield stress identified 
in rheological flow curves. Furthermore, we suggest that an alternating shear of frequency 
CO and amplitude 7 should lead to a dynamic phase transition line in the (a),y) plane, from 
a 'flowing' to a 'non-flowing' phase. Finally, we apply our results to crystals, show the 
corresponding microscopic process leading to flow and discuss possible experimental inves- 
tigations. 



We clearly know how to make the difference between a solid and a liquid in everyday life. 
However, the question of whether solid and liquid are distinct states of matter is far from 
being conclusively answered in the literature. A difficulty in addressing this issue is that even 
ordinary liquids may look solid, as one learns when diving into the water from an height 
of, say, 10 meters without being a pro. Actually, supercooled liquids and soft matter, such 
as emulsions, foams, and pastes behave like bona-fide solids, if probed on timescales that 
range from microseconds to hours. This phenomenon, called visco-elasticity, is a property 
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of most liquids: they respond as elastic solids on short timescales and as a real liquid on long 
timescales (for an amusing illustration see the pitch experiment UJ). Of course the meaning 
of long and short depends on the system at hand: the viscous-elastic crossover takes place 
on timescales of the order of 10^'^ s for water[^and on the order of hours for supercooled 
glycerol close to 190K. The question we address in this work is whether all known solids are 
actually like that, i.e. liquid-like on large enough timescales. This is the conjecture of some 
rheologists, as nicely argued in [2|. On the other hand, general theoretical results, see e.g. 
t3i , argue that rigidity is a consequence of breaking translation invariance (in all directions): 
states of matter characterized by this spontaneous symmetry breaking, like crystals, must be 
rigid and have a finite elastic shear modulus. 

How can we reconcile these two apparently antitethic claims? First, it is important to 
remark that solids in nature contain often a large quantity of non-equilibrium defects that 
allow them to flow rather easily. The main example is provided by dislocations in crys- 
tals. They cost infinite energy]^ and, hence, they are absent in three-dimensional crystals at 
thermal equilibrium. Although equilibrium crystals can be prepared ll5]|6l, typical growth 
procedures lead to imperfect crystals characterized by a dislocation density that ranges from 
10*^-10' cm . When a shear stress a is applied, dislocations move with a velocity pro- 
portional to the stress and lead to what is called plastic flow 17 |. The experimental examples 
studied in |2] are probably all characterized by some kind of non-equilibrium defects. On 
the other hand, the theoretical results cited in |31 refer to equilibrium states of matter that 
break completely translation invariance, which we will call in the following "perfect solids". 
One possibility would be therefore that perfect solids are rigid but that non-equilibrium ones 
flow. Instead, here we show that even perfect solids do flow on long enough timescales if an 
infinitesimal macroscopic shear stress is applied. Despite this property, we also show that 
perfect solids are different from ordinary liquids because their viscosity diverges for vanish- 
ing shear stress with an essential singularity — contrary to the finite viscosity characteristic 
of ordinary liquids. The role of the thermodynamic limit is a subtle and important issue. In 
fact, it is clear that any finite system cannot be rigid and will flow if one waits a long enough 
time, see e.g. ID. The crucial question is whether or not this time scales with the system 
size. In the former case real solids would flow but this would happen on timescales too large 
to be relevant for any practical purposes, maybe even larger than the age of the universe. We 
shall show that instead the timescales for flow are long but finite in the thermodynamic limit 
and attainable in experiments. 

Our starting point is similar to the usual nucleation arguments for metastable phases. 
Consider a solid of size (where d is the spatial dimension), — it could be an equilibrium 
phase such as a crystal, a quasi-crystal or an ideal glass, if such a thing exists; or even a 
metastable state such as a crystal with many entangled dislocations. We apply a very small 
but finite shear stress ct, see Fig. [T] After a very rapid equilibration time the system will 
deform by an amount y= o/Gp where Gp is what is called in practical applications the 
elastic shear modulus. This deformation leads to an increase in the free energy of the system 
that according to elasticity theory is given by: U' This implies that the deformed solid 
has a free energy extensively higher than the one of the equilibrium undeformed state and, 
as a result, that the deformed solid is metastable and has a finite life-time, except if the free 



' Therefore, in the diving example water 'looks' sohd but not because of visco-elasticity. See fsl for a 
detailed explanation of the physics of splash. 

^ The energy of a single dislocation diverges so strongly in the themiodynamic limit for dimension higher 
than 2 that the equilibrium dislocation density is zero fj]- 
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Fig. 1 (Color online) Schematic representation of the deformation due to a shear stress. The amount of 
deformation is measured by the angle 7 which is proportional to the shear stress a. The proportionality 
constant is the inverse of the elastic shear modulus G. 

energy hairier to escape from this metastable equilibrium is infinite - which is not the case, 
as we shall show in the following. 

The introduction of a 'static' modulus begs an explanation, since we shall argue in what 
follows that even perfect solids flow under finite stress. We shall show that under the effect 
of a stress, the system responds in two time-steps: a rapid one corresponding to the elastic 
response inside the metastable deformed state, followed by a slow flow on large timescales. 
The value Gp{t) we use here is defined as the elastic shear modulus obtained at times long 
compared to the elastic deformation, but short compared with the stress-dependent nucle- 
ation time|2] The two-timescale response of course applies (by definition) to viscoelastic 
systems, but we shall argue here that it holds also for true solids at long enough timescales: 
the difference with viscoelasticity being that for a solid the long-time flow rate vanishes 
faster than a at small stresses, thus implying an infinite viscosity in the limit of zero stress. 

Let us now introduce our main argument by focusing on a cavity of radius R inside 
the deformed solid. Imagine reshuffling all particles inside the cavity from the deformed 
to the undeformed configuration. For example, in the crystal case this is shown in Fig. [2] 
Then, while still maintaining the applied shear stress a, we let the solid evolve dynamically 
from this configuration. After a short equilibration time the system reaches a metastable 
equilibrium. The free energy difference between this state and that without cavity reads: 



where the first term corresponds to the free energy cost one has to pay due to the mismatch 
at the surface of the cavity between deformed and undeformed state. Q (a) is like a surface 
tension and can in principle depend on a. For small a, the limit we are interested in, the 
dependence on G matters only if 12 (ct) does not go to a constant and vanishes]^ We will see 
that this is indeed what happens for crystals. 

AF{R) first increases as a function of R, has a maximum at R* = Gpnja) 
then decreases for R > R*. This implies that the nucleation of a droplet of the undeformed 

' This is in strict analogy with the susceptibility in a ferromagnet: a small field applied against the magne- 
tization will induce a linear response at intermediate times, but will eventually destabilize the phase. 

** We consider particles interacting with a short range non singular potential, therefore it is sure that £2{a) 
has to be bound by a finite positive constant for C7 — > 0. 



AF{R) =R''-^n{a)-R' 




(1) 



Fig. 2 Defects structure around a relaxed cavity embedded in a slieared system. Tlie gedanken experiment 
discussed in the main text consists in reshuffling all particles inside the cavity from the deformed to the un- 
deformed configuration. Here we show the case of a two-dimensionnal hexagonal crystal for a strain 7 = 0.4. 
The interaction between particles is of hard core type with a radius O.35o, where a is the lattice spacing. 
Particles are plotted as disks with radius a/2. White (empty) particles have six neighbours and thus satisfy 
the hexagonal order, whereas red and blue ones correspond respectively to particles having seven and five 
neighbours, which are elementary topological defects called disclinations. Note that a 5-7 dipole of disclina- 
tions corresponds to a dislocation, the Burgers vector of which is perpendicular to the dipole axis. Here, the 
cavity has a radius bigger than R* , the critical radius of a relaxed nucleus in a such system, see text. 



solid of size R is unfavorable for radii less than the critical value R* . In contrast nuclei 
larger than the critical one of size R* will grow since they are free-energetically favored 
Thus, the free-energy barrier for nucleation via a spherical droplet is given by AF{R*) = 
2Cp(d-\) ' ^^'^ '■'^'^ clearly an upper bound. It is necessarily finite for any finite a, 
although it diverges for ct — )• 0, as is usual in metastability problems (here ct plays the 
role of the external field inducing metastability). Depending on the value of a and system 
parameters, either different nuclei are already present in the system with a density p„ ~ 
1/ {R*Y exp{—AF{R*)/kBT) and grow with a finite speed or there are none and one has to 
wait a time of the order Tq exp {A F {R* )/kBT) to nucleate one that then will grow and invade 
the system. In both cases, for very small a, the characteristic timescale for destabilizing the 
metastable deformed solid is given by T* = Tq exp{AF {R* } / ksT) where Tq is a microscopic 
timescale that will be discussed later. 

As a consequence, after a time of the order T*, the deformed solid goes back to the 
undeformed state. If, however, the external shear stress ct is continuously applied, the solid 
deforms again and the entire process repeats. Although a more refined analysis of the flow 
process is needed (see below for a first attempt), these simple arguments already suggest 
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that the solid deforms with a shear rate 



7 CT a 

— = = exp 

T* GpT* GpTQ 



n \ 2d-2 \ 



where cJv(o') = Gp yi^^^T^ Wrj effective scale for the stress and c is a con- 

stant, which equals k/2> in three dimensions. If Q[a) does not vanish for cr —>■ nei- 
ther does Oy', if instead Q{a) vanishes, then the previous law connecting the shear rate 
7 to the shear stress CT changes and the exponent 2J — 2 is renormalized. In any case, as 
long as CTy vanishes less rapidly than linearly in the stress (for crystals we will show that 
Oy CT^/^(—log(CT/G/)) the previous equation leads to an essential singularity for the 
viscosity of the system: 

CT / / CT \ 2d—2\ 

,^ = _=G^Toexp(^c(^) j (3) 

This is the main result of this work: a perfect solid does flow but with a viscosity that di- 
verges extremely fast when CT — >^ 0. The conclusion of the previous arguments is that what we 
call in practical applications "elastic shear modulus", i.e. Gp, is, as described above, only an 
effective timescale-dependent quantity corresponding to the value of the very long (but not 
extending until zero frequency) plateau displayed by G{(o) I15II . Similarly to visco-elastic 
liquids, perfect "solids" behave elastically on timescales much less than T* and eventually 
flow on timescales larger than T*. However, the flow properties for small shear stress are 
very different: instead of the singular relation (|2]i, ordinary liquids are characterized by a 
linear Newtonian relation leading to a finite viscosity at CT = 0. 

Let us note that the metastability of the stressed state is closely analogous to the metasta- 
bility of current-carrying states in superconductors and superfluids 1 10], and of loaded states 
in solids IllJ. Indeed, the nucleation argument applies to all situations where a system breaks 
a continuous symmetry, and acquires a rigidity with respect to long-wavelength deforma- 
tions along this symmetry f3l: whenever the physical situation is such that the "twist" along 
the zero-modes per unit length is finite in the broken symmetry state, the phase becomes 
metastable. Such is the case for superfluids, superconductors, solids and smectics. In the 
context of two-dimensional solids our argument is very similar to the one already devel- 
oped in I12II13I in order to obtain estimates of the viscosity in the strong strain regime. Our 
purpose is, however, quite different: we give a weaker argument that just gives us an upper 
(finite) bound for the viscosity, but such that it may be applied in any dimension, and without 
invoking the nature of the order responsible for solidity. This allows us to include (putative) 
ideal glasses whose spatial order is poorly understood at present. It is reassuring that in the 
case of two dimensional solids, where a detailed analysis can be performed 112] , our upper 
bound is just slightly sub-optimal|^ 

At this stage, it is important to address an apparent paradox: there are rigorous expres- 
sions relating the elastic shear modulus to correlation functions evaluated in the equilibrium 
(unstressed) solid I16II17I . These are used in simulations to obtain quantitative values of the 
elastic moduli and were shown to be different from zero. There is actually no contradiction 
with our results, but the reason is subtle and worth discussing in detail. These equilibrium 



^ Our argument gives log T « — log y <x l^°s{°/Op)] jj^g limit of small shear stress, whereas the true 
behavior is log7°<: [\og{a/Gp)] and corresponds to an unbinding mechanism of the dipoles of dislocations 
that are present in the system in the low temperature crystalline phase. 
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expressions of the shear modulus are obtained via linear response and a fluctuation dissi- 
pation relation expressing the response of a crystal to a shear deformation in terms of the 
correlation functions of the undeformed crystal. This linear response is justified only to the 
extent that the stress does not induce a macroscopic change, and this will happen either if 
strain is microscopic (y not of order one), or if it acts for a time T < T* only. Either cases 
yield a non-zero modulus Gp within the deformed metastable state. As for the general the- 
oretical results, cited in the introduction |3 |, stating that the existence of a shear modulus is 
a consequence of breaking translation invariance in all directions, a careful analysis shows 
that the breaking of translational invariance implies the finiteness of Gp (defined for times 
such that the system stays within the metastable non-flowing state), and not of the true shear 
modulus associated with arbitrarily long times and finite stress. Spontaneous thermal fluc- 
tuations will not make a true solid flow - hence the broken translational invariance, but an 
external strain of the order of the system size will. Again, this is in strict analogy with a 
ferromagnet, where spontaneous magnetization fluctuations will not destabilize a phase, but 
a finite external magnetic field will. 

The argument above not only implies that nucleation will make an equilibrium solid 
under finite stress start flowing, but also that it will preclude it from stopping. In the case 
of a crystal, as we shall discuss below, the flow is facilitated initially by dislocations that 
are nucleated. As is well known, the dislocations evolve and may entangle, and their motion 
become restrained. However, the out of equilibrium phase of interacting dislocations cannot 
sustain a finite stress forever, since the nucleation of an unstressed phase is also possible in 
that case. What one may observe in such a situation is a first regime in which the defects 
proliferate and entangle, and a second regime, with much slower flow, in which the system 
responds in a manner more similar to that of an amorphous solid. Needless to say, this also 
applies to a situation in which there where some nonequilibrium defects present in the initial 
configuration. 

We now apply the previous general arguments to the initial flow of perfect crystals. In 
order to find the a dependence of 12 (ct) we consider a cylindrical cavity of height and radius 
equal to / = ^ = and whose symmetry axis is perpendicular to the shear direction (a is 
the lattice spacing). In this case, our construction boils down to creating a single dislocation 
loop with radius a/y and Burgers vector of length a. Both the loop and the Burgers vector 
are parallel to the basis of the cylinder. Creating the dislocation loop has an energy cost 
4^(i'-v) [^°8(' /'') where A is a constant of the order of one and v is the Poisson ratio (T). 

However, stress is also released, which leads to an energy gain proportional to = ^^fi 
171 . Thus, on the scale Z = ^, we have found a microscopic expression for the energy gain 
and energy loss terms introduced previously. By equating the generic expression Q.{a)f' to 
its value for Z = ^ we find Q{a) = 4;r(i'"v) [~ log(<^ /Gf ) — A] where c is a constant, whose 
value depends on the geometry of the cavity and which we will take equal to one in the 
following. Assuming v :^ 1/3, which is a typical value for solids, one obtains the following 
laws for crystals in the limit of very small a: 

T* = Toe* , 77=GpTo/ <^ , with a^. = -P—, (4) 

where is a constant whose value is close to one and is system-dependent; Gy is an effective 
stress scale. Note that since generically Gp :^ kgT /a^ |3|, the typical value of Oy is of the 
order of Gp (or slightly larger). The value of K* for crystals turns out to be of the order of 
— alog{(J /Gp)/ {(J /Gp). Thus, roughly, the critical nucleus is formed by a few dislocation 
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loops, whose nucleation has already been discussed in the literature, see e.g. fTl llSI . This 
also allows one to estimate the microscopic timescale. Whether the critical nuclei are already 
present in the system or one has to be nucleated, Tq should be roughly given by the time to 
nucleate a kink-antikink pair in a dislocation line (in 3D), which is of the order of the Debye 
time multiplied by exp{2W / ksT) where W is the kink energy. This is in fact the typical time 
to nucleate a microscopic dislocation loop and also the time on which a dislocation moves 
of one lattice spacing. 

In the previous arguments, the non-equilibrium steady flowing state was described in a 
very phenomenological way. In the following we investigate in more detail the flow process 
and endeavor to obtain some predictions for the flow curve of a "solid". For not too small 
shear stresses, there is a finite density of critical nuclei]^ and one has to take into account 
their growth and their continuous deformation due to a constant shear rate. Forgetting the 
continuous deformation for the time being, one finds that the growth is a two-stage pro- 
cess: critical nuclei first expand driven by a free-energy gain, until they start to overlap. At 
this point the growth becomes a coarsening-like process |19|. Since this is expected to be 
much slower, we will neglect the latter in the following. We model the system as divided in 
mesoscopic blocks characterized by a certain value of the stress a. The block sizes are of 
the order of the distance between critical nuclei. We assume that the stress inside the block 
increases with time due to the constant shear rate. A critical nucleus can appear inside the 
block with a rate I /t* (cr) (see eq. ([2|). This leads to a subsequent decrease in the stress that 
we model with a transition a — )■ OP] The equation for the evolution of the distribution of 
block stresses, P{a,t), is therefore: 

5,P(a,0 = -Gp75^P((j,0-^P(a,r) + 5(a)^^°°^;(j'^J^P((j',f) (5) 

The stationary solution of equation ^ can be computed exactly. It reads P(ct) = Cexp ^— Jq 

with C a normalization constant. From P{o), one can obtain the mean value of the stress in 
the system and the corresponding flow curve, see Figure [3] Note that in the limit of small 
y, we find back the relation l[2| since the integral over P(cj) is dominated by a's such that 
g^y^t^^j °^ 0{V). For high shear rates, when the typical distance between nuclei is rather 
small, a rearrangement of the stress inside one block will change somehow the stresses in 
the nearby blocks. We have taken this effect into account by modifying the previous equa- 
tion following the approach developed in |20,21|. The results will be published elsewhere 
|22J . They lead to a refined description of the flow curve but do not alter our conclusions. 
Note that for very small shear stresses, the flow is due to the nucleation of just one critical 
droplet of unstressed material and its subsequent growth until reaching the whole system. 
In this regime the coarsening-like regime will not be present and the previous description 
is expected to apply even better. One important conclusion of the analysis of our simple 
model is that there is no yield stress since the system always flows. However, the ultra-slow 
logarithmic decay of the stress a <^ iog{Gf) easily fool the eye and lead to the definition 
of an effective yield stress when only a limited range of y are available, see Fig. [3] Thus, 
our reasoning has a consequence for the Jamming Diagram of Liu and Nagel |23], since we 
have seen that the truly non-flowing phase collapses to three subsets of the 0=0, the T = 0, 
and the P = oo planes, respectively. 

* From a purely statistical mechanics point of view, if one considers an infinite system then there is always 
a finite density of critical nuclei. 

^ Clearly this is not an instantaneous process and a more refined model should account for the finite time 
it takes to the nucleus to growth and invade the block. 
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Fig. 3 Flow curve, C7 as a function of / resulting from equation |5|. The stress fj decreases extremely slowly 
(see text) when y — > 0: it can be interpreted in terms of an apparent yield stress in experiments and simulations. 

A very intriguing question arises when one considers oscillatory shear with frequency 
03. Ferromagnetic systems subjected to an alternating magnetic field have a dynamic phase 
transition from a low-frequency regime in which nucleation destroys order, to a high fre- 
quency regime with spontaneous average magnetisation 1.24^ 14J . It seems likely that an 
analogous dynamic phase transition, characterized by the vanishing of the amplitude of the 
Bragg diffraction peaks, occurs for the flow of solids. One then expects that the density 
profile of a solid will perform periodic motion, and would hence stroboscopically stay put, 
under periodic shearing throughout a phase of low amplitude and high frequency. On the 
contrary, for frequencies less than 1/T*(y), the nucleation mechanism we discussed above 
will become effective and will lead to a translationally invariant steady state. This also ap- 
plies to other kinds of solids, such as quasicrystals and ideal glasses (if they exist). Within 
the flowing phase, the stroboscopic density, averaged over many cycles, will be spatially flat, 
while it will have some modulation (quasiperiodic, amorphous) in the non-flowing regime. 

Our results unveil a general possible mechanism inducing flow of solids (e.g. crys- 
tals, quasi-crystals or glasses). Thus they provide upper bounds on the relevant timescales. 
Whether or not there exist more complicated and faster process is an open questioij^ In any 
case, it is interesting to discuss the possibility of observing in experiments the flow phenom- 
ena discussed in this work, in particular for three dimensional crystals. Using the value of the 
kink energy estimated for example for Aluminium, W ~ O.QlSeV I25L we find Tq °<: IQ^^s at 
ambient temperature (working at higher temperatures would diminish this value). Using our 
previous results, in particular Fig. pi we find that imposing a stress of the order of lO^'Gp 
lead to deformations of the order of one on timescales of the order of 10^.?. This is clearly a 



For example, for crystals, depending on the way the solid is deformed, the Nabarro-HeiTing creep mech- 
anism could play an important role. 
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regime that can be investigated in experiments. Usually, real crystals undergo fracture with 
this kind of values of o, but that is because they contain flaws and non-equilibrium defects. 
The pre-requisite to discover the effect predicted in this work is therefore to prepare perfect 
crystals. This is not an easy procedure but it can be done |5 , 6]. Colloids are other interesting 
systems to test our results. Micro-rheological experiments in which the motion of particles 
is visualized by confocal microscopy have been already performed. Actually, it has been 
possible to study and visualize nucleation of a single dislocation loop 1 18|. This strongly 
suggests that observing nucleation of the unstressed state and a dynamic phase transition 
under alternate stress is within reach. The microscopic timescale for colloids is of the order 
of ?ns; as a consequence, Tq becomes much larger and this pushes all timescales roughly 9 
orders of magnitude up. Thus to obtain yoc 10^^, one should impose a stress of the order 
of 0.3 Gp. This drawback is compensated by two advantages: for hard-sphere colloids, one 
can rather easily prepare perfect crystals |18| and even imposing deformations 70c 10^' 
should not lead to fracture. Clearly simulations can also be a very useful tool to unveil the 
existence of nucleation of unstressed regions inside a stressed solid. Although nucleation is 
a rare event and, therefore, difficult to obtain in a finite size/time numerical simulation, one 
could already study the growth of the unstressed regions by starting from an initial condition 
with the nucleus present. 

Finally, the generalisation of our study to glassy liquids is worth studying, see [261 for 
a work in this direction. Close to their glass transition they are to a certain extent half-way 
between solids and liquid. The repeated nucleation of unstressed regions could be relevant 
to explain rheological properties for small applied shear stress close to the glass transition, 
where an apparent yield stress emerges. 
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